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ABSTRACT 

We study aspects of confinement in two deformed versions of the AdS/CFT 
correspondence - the GPPZ dual of M = 1* Yang Mills, and the Yang Mills* Af = 
dual. Both geometries describe discrete glueball spectra which we calculate numer- 
ically. The results agree at the 10% level with previous AdS/CFT computations 
in the Klebanov Strassler background and AdS Schwarzchild respectively. We also 
calculate the spectra of bound states of the massive fermions in these geometries and 
show that they are light, so not decoupled from the dynamics. We then study the 
behaviour of Wilson loops in the lOd lifts of these geometries. We find a transition 
from AdS-like strings in the UV to strings that interact with the unknown physics 
of the central singularity of the space in the IR. 
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1 Introduction 



The AdS/CFT correspondence [TJ E] has led to the development of dual string theory (su- 
pergravity) descriptions of a variety of large N gauge theories. As originally stated the duality 
applied to gauge theories with maximal supersymmetry including, in the case of four dimensions, 
M = 4 Super Yang Mills. An interesting task is then to try to move in the direction of gauge 
theories which are more phenomenologically applicable. This necessarily means looking at cases 
with less than the maximal supersymmetry of the original duality. Several examples of such 
dual geometries have been constructed including the non-supersymmetric AdS Schwarzchild ge- 
ometry [5 , and the M = 1 geometries of Klebanov Strassler j3] and Maldacena Nunez jH] . It is 
now well understood how to compute the glueball bound state predictions [7J and Wilson loop 
behaviour [Hj in such theories and these geometries show confining behaviour. 

In this paper we will consider some of the simplest interesting deformations of AdS, which 
will provide results for glueballs and confining strings that test the systematic errors within such 
AdS approaches. The geometries we will study are asymptotically AdS with deformations in the 
interior that can be interpreted as including relevant (mass) operators in the gauge theory. Thus 
the theories are strongly coupled M = 4 super Yang Mills in the UV but have less symmetry in 
the IR. The simplest method for constructing such deformations ^0] is to use five dimensional 
M = 8 gauged supergravity where the scalar fields act as source and vev for operators in the 
field theory. Lifting the solutions back to lOd ^1] is somewhat involved but the lifts of the 
two geometries we will study are in the literature [121 HHj (in fact the M = 1* lift is not quite 
complete). 

The first of the geometries we study is the GPPZ M = 1* solution where one of the 
scalars of Af = 8 supergravity has a non-trivial profile corresponding to a mass term for three 
of the four fermionic fields of M = 4 super Yang Mills. The solution is supersymmetric so 
the 6 scalars are also massive. The IR of the theory has just a gauge field and a gaugino. 
Glueballs in such geometries correspond to excitations of the dilaton field in the background. 
The potential for such fluctuations is a well [T3] and therefore gives rise to discrete solutions, 
showing that the theory has a mass gap and is confining - here we numerically compute the 
solutions and find the mass spectrum. The results match to two significant figures the predictions 
[Ho] from the Klebanov Strassler [S] construction, which should describe a similar IR gauge 
theory. This encouragingly suggests that the UV completions of the theories, which are different, 
are unimportant to the glueball spectrum. In the GPPZ solutions one can also switch on a second 
scalar corresponding to a real gaugino condensate [H]. Previous analysis has suggested that 
the solutions with too large a condensate are unphysical. From the glueball analysis we confirm 
the existence of an upper limit on the gaugino condensate parameter since above it the glueball 
potential becomes unbounded from below. Within the physical solutions the value of the gaugino 
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condensate only changes the glueball masses by of order 10%. It is also straightforward, following 
the analysis of [T7j . to calculate the masses of the bound states of (AA) gauginos. Since there 
are a set of solutions with different gaugino condensate the gravity dual describes a moduli 
space and therefore we find a massless bound state plus excitations. The degeneracy of the 
vacua is a symptom of large N and is unfortunately unhelpful for comparison to finite N lattice 
computations ^H] where there is no massless state. One might hope to compare our glueball 
results to future lattice results though. Finally we test the decoupling of the UV degrees of 
freedom by calculating the masses of the bound states of heavy fermions; we find they are of 
comparable size to the glueball masses so the decoupling is far from complete. This makes the 
insensitivity of the glueball masses to the UV physics all the more remarkable. In principle one 
should take into account mixing between these states and those we consider. Such an analysis 
can be found for the M = 1* theory in [22] ■ Comparison to those results shows that the mixing 
again only effects the glueball masse at the 10% level. 

Our second geometry is the Yang Mills* dual JSI which in five dimensions has a single scalar 
switched on, corresponding to an equal mass for all four fermion fields of M = 4 super Yang 
Mills. Brane probing the lOd lift ^3] shows that the six scalar fields also have a mass and 
therefore the IR of the theory is just a non-supersymmetric gauge field. We find the potential 
relevant for the glueball spectrum of the theory, show that it is a bounded well and compute 
the discrete glueball spectrum. The spectrum matches at the 10% level with previous estimates 
of the spectrum from AdS Schwarzchild [Zj. This is again a surprising level of insensitivity to 
the UV completion since the AdS Schwarzchild geometry describes a 5d gauge theory in the 
UV whilst YM* returns to 4d M = 4 super Yang Mills. We also calculate the fermionic bound 
state masses and find they lie close to the glueball spectrum showing the limitations of the 
decoupling. 

Although these deformed geometries are easy to construct and work with we must be careful 
because they are singular in the deep interior. Naively one would conclude that there is some- 
thing badly amiss! However, in a similar construction of the M = 2* theory ^01 1201 EH 1221 12H| 
the interior is also divergent. There brane probing has revealed that the singularity matches 
to the pole of the one-loop coupling of the gauge theory (the enhangon mechanism [21] )• This 
is a possible interpretation of the singularities here too. It is also possible that the configura- 
tions are analogous to the Polchinski Strassler M = 1* construction ^B] in which the core of 
the geometry contains some complicated fuzzy expansion of the D3 branes although identifying 
such an object from the geometry appears hard. Further one might expect that the strong 
coupling might trigger additional operators to switch on and these might act to smooth the 
singularity. This is the mechanism that leads to smooth geometries in the Klebanov Strassler 
and Maldacena Nunez geometries. Thus our geometries may well be incomplete. Nevertheless 
they do have the correct UV behaviour and conformal symmetry breaking. We might hope that 
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typical strong coupling phenomena will be triggered by a large but finite gauge coupling and 
that the geometry very close to the singularity will therefore not play a role in the geometry's 
predictions. The glueball computations support this conclusion since the bounded well that 
controls their mass forms before the singularity is reached. 

To investigate this further we also study Wilson loops in these geometries. The Yang Mills* 
lOd lift is the most complete so we concentrate on that case. A probe fundamental string in 
the geometry is associated with the interaction between a test quark anti-quark pair in the 
background gauge configuration jH] . In pure AdS the wider the quarks are separated the deeper 
the string penetrates into AdS. The action of the string looks like a constant (corresponding to 
the quark mass) plus a term determined by conformal symmetry that goes as —1/L where L 
is the quark separation. In the Yang Mills* geometry the solutions of the equations of motion 
of the string behave very differently. Asymptotically the geometry is AdS so for quarks that 
are close the behaviour is the same. However, as the string samples further into the deformed 
space the Euler Lagrange equation solutions describe a second solution in which the quarks are 
closer together again. We show that this solution has higher action than the AdS-like solution. 
Naively there appears to be a maximum length string (we suggested this interpretation in the 
first version of this manuscript). To understand the solutions better we study the action of a set 
of curves linking two quarks at each fixed separation. In fact a simple sine wave configuration 
linking the quarks is sufficient to reveal the role of the different solutions, even though this set of 
curves only contains an approximation to the true solution of the Euler Lagrange equations. We 
find that the second solution is in fact a global maximum of the action. Strings that penetrate 
deeper into the core than that configuration have falling action and the string in fact wants to 
fall into the singularity. We find a phase transition as the quarks are separated where AdS-like 
strings are replaced by strings entering the core. This is unfortunate since it means the physics 
described by the geometry can not be completely described by just the supergravity. A stringy 
resolution of the singularity is needed. One can imagine that all might yet be well since if 
there is some fuzzy brane construction then the set up might be consistent with confinement 
by mechanisms similar to those suggested in the Polchinski Strassler M = 1* set up ^H]- Our 
knowledge is clearly incomplete though and this should be taken into account when considering 
the glueball spectrum results. 

We begin in the next section with a review of the formalism for constructing deformed AdS 
theories and looking at fluctuations about such geometries. In section 3 we will study the 
glueballs of M = 1*, and in section 4 the glueballs of Yang Mills*. In section 5 we turn to the 
analysis of Wilson loops in these geometries. 
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2 Deformations and Fluctuations 



We briefly review the formalism for constructing deformed AdS geometries in 5d supergravity 
and for studying the fluctuations about such backgrounds. 

2.1 Deforming AdS 

The AdS/CFT correspondence PH E| maps supergravity fields to operators in the field theory. 
The theories we consider here correspond to deformations of M = 4 Super Yang Mills by mass 
terms for the four fermions of the theory. These operators correspond to scalar fields of the 
5d supergravity theory. The dependence of the supergravity field on the radial AdS coordinate 
encodes information about the RG flow of the field theory operator. 

Thus we look for solutions involving only scalar fields and allow the scalars to vary only in 
the radial direction of AdSs The other non trivial field is the metric 

ds 2 = e^^Wttep + dr 2 , (1) 

where /i = 0, . . . , 3, r is the radial direction in AdSs, and A(r) — > r as r — > oo so that the space 
becomes AdSs at infinity. With an appropriate parametrization the scalar lagrangian can be 
written as C = |(<9A) 2 — V(X) and the equations of motion for A and A are [TO] 

X"+AA'x' = ^, 6A' 2 = \' 2 -2V. (2) 

CM 



For large r, where the solution will return to AdSs a t first order and A — ► and V — > ^rj-A 2 , 
only the first equation survives with solution 

A = ae-^ r + be~ Ar (3) 

a and b are constants, while the conformal dimension A of the field theory operator is related 
to the mass of the supergravity scalar by 

M 2 = A(A-4). (4) 

a is interpreted as the source for a field theory operator and b as the vev of that operator, since 
e r has conformal dimension one. 

If the solution retains some supersymmetry then the potential can be written in terms of a 
superpotential 



dW 



--\W\ 2 (5) 



3' 



and the second order equations reduce to first order 



A = ~-=r, A = — W. 6 
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2.2 Linearized Fluctuations 

The supergravity scalars characterised by a non zero 6 in (J3J) map to field theory operator 
expectation values. We can therefore look for linearized fluctuations of a given scalar in a 
background, corresponding to bound states in the field theory associated with that operator. 
We follow the analysis in |H Q E] • We look for normalisable solutions of the linearized wave 
equation 

d,>W=gsrW* = ^5$, (7) 
that behave as plane waves at the boundary 

5$ = ^j( r )e~ ikx . (8) 

To determine the spectrum of the fluctuations, a standard procedure in AdS is to reduce the 
above equation to a solution of a Schroedinger problem [U El HI] • Making the change of coor- 
dinates (r — > z) 

and rescaling 

$ _> e - 3A/2 v, (10) 

equation (JZj) takes the Schroedinger form 

(-^ + f/( 2 ))#) = M 2 #), (11) 

where 

17 = - A" + -(A') 2 + e 2A S^- (12) 
Primes now denote differentiation with respect to z. If this potential takes the form of a 

bounded well then there are discrete solutions (with k 2 = —M 2 ) for the linearized fluctuation 

indicating confinement of the fields into discrete bound states. 

Since we will work throughout in the z coordinates, we re-express the results of the previous 

subsection in these new coordinates. The second order equations of motion for the background 

deformation become 

A" + 3X'A' = e 2A ^, 6A' 2 = \' 2 - 2e 2A V. (13) 
In the UV {z — > 0) the behaviour is 

z = -e~ r , e 2A = \, X = a(-z) 4 - A + b(-z) A , (14) 

z A 

where, again, a is interpreted as a source for an operator and h as the vev of that operator since 
z has conformal dimension -1. Finally the first order equations © now read 



5 



3 The AT =1* Geometry 



The M = 1* theory is M = 4 super Yang Mills with equal mass terms for the three adjoint chiral 
superfields leaving just the vector multiplet massless. The UV of the theory is conformal whilst 
the IR behaves like M = 1 Yang Mills generating a gaugino condensate dynamically in the weak 
coupling limit. At large N, where we will be working, the UV is also strongly coupled so the 
chiral multiplets can not be considered fully decoupled. The theory at large N has been studied 
in and it has been shown to have a set of discrete vacua differentiated by the magnitude of 
the gaugino condensate which is real. 

The gravity dual of this theory was found by GPPZ . In the 5d supergravity theory one 
must identify the scalars dual to the the dimension 3 operators 

4 

m = EV#i, O ir = M 1 . (16) 

i=2 

These corresponds to two scalars m and a in the 10 of the SO (6) gauge symmetry of the 
supergravity theory (the global SU(4) fi symmetry of the gauge theory). 

The GPPZ flow tuned by these two scalars preserves an SO (3) flavour symmetry. The 
supergravity fields that are singlet under this residual 5*0(3) can be organized in one gravity 
multiplet and two other hypermultiplets 

The potential for a and m is 



i /in y in 

V = -^[cosh 2 (^) + 4cosh(^)cosh(2a) - cosh 2 (2a) +4], (17) 



8 L v v^ 7 K V3 
and can be obtained from a superpotential 



2 ui 

cosh(^^) + cosh(2<r) 



(18) 



V3' 

as expected since the solutions should maintain M = 1 supersymmetry. We can then work with 
the first order equations 

% = -^sinh(2a), (19) 
dm v3 a , ,2m. , , 

& = -- esmh( T3 h (20) 

^ = ie' 4 [cosh(?|)+cosh(2a)]. (21) 

Equations (|19H21j) can be solved analytically jH] , and the result is a family of solutions parametrized 
by the boundary values of the fields as z — » (UV) 

m — >• —az, a^—bz 3 , A—>—\og\z\. (22) 
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From (|14|) . m corresponds to a mass for three of the adjoint fermions and a is the gaugino 
condensate. In this paper we do not need the analytic expression for the flows, so we will solve 
the first order equations numerically subject to the boundary conditions (j22J- 

Solutions exist for all initial choices of a and b but the resulting geometries are singular. 
In [T2j it was argued that only a subset of these flows should be considered physical. There it 
was required that the supergravity potential evaluated on the solution should be bounded from 
above. This condition restricts us to flows with initial conditions 

b < 3~ 3/2 a (23) 
» 0.19a 

The gravity dual therefore predicts a moduli space of vacua with varying real gaugino condensate 
up to some maximum value (0.19m) which is plausibly the iV — > oo limit of the field theory 
result discussed above. 

We will see further evidence for this bound on a shortly. For convenience in what follows we 
will set a = 1 - that is we set the scale of all mass operators in terms of the chiral multiplets' 
mass. 

3.1 Glueballs 

The mass spectrum of bound states can be obtained from fluctuations of the dual scalar fields 
by finding the eigenvalues of ([lip. We first calculate the mass spectrum of glueballs with 
quantum numbers J PC = ++ associated with the correlators (Tr F 2 (x)Tr F 2 (y)). The operator 
O = Tr F 2 is dual to the dilaton, so we need the Schroedinger potential ()12)) in this case. 
The dilaton does not contribute to the supergravity scalar potential, so that the Schroedinger 
potential for the glueballs is just 

U= 3 -A" + l(A') 2 . (24) 

Actually the dynamics of the eight scalars of the full supergravity background is quite compli- 
cate and the dilaton, though being inert (that is having no radial dependence in the background 
solution) couples to the other scalars. To find the complete Schroedinger-like equation for the 
dilaton one has first to diagonalize the equation of motions for the system of scalars, as in , 
where an analytic solution has been found for all the scalars. Bound state spectra for some 
scalars had been previously calculated also in (2E1 123 I2E 123 E2j • 

In order to obtain a discrete spectrum with a mass gap U must be bounded below. We plot 
the potential for varying b in Figure 1. The potential is indeed bounded for all physical flows 
satisfying the bound ()23|) . The unphysical flows with b > 1.9a, which violate (|2*3*j) . have an 
unbounded glueball potential. Thus we have further evidence for the bound. 
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The eigenvalues of (JTTJ), and hence the glueball mass spectrum, can be easily obtained using 
the numerical shooting method. Since we want the solution to match to the operator Tr F 2 we 
set the UV {z — > 0) boundary conditions on ip to be 

m = j v (25) 

and then numerically solve for various M seeking regular solutions in the IR. 

The results for 6 = and b = 0.19 are shown in Table 1. The ++ mass is not a prediction, 
but just sets the strong coupling scale, so we normalize the lightest glueball state to 1. The 
gaugino condensate's value has little influence on the glueball masses (at most of order 10%). 

We also quote the glueball mass spectrum obtained from the M = 1 Klebanov Strassler 
model ESI- The Klebanov Strassler dual also describes Af = 1 Yang Mills in the IR but in 
the UV there are extra flavours that participate in a cascade of Seiberg dualities. The glueball 
predictions agree very well with the M = 1* results, suggesting the UV completion of the theory 
is relatively unimportant to the glueball spectrum. The Maldacena Nunez [Hj dual also describes 
M = 1 Yang Mills in the IR and the glueball spectrum has been recently analyzed in [HB] • That 
dual though becomes strongly coupled in the UV so the computation is less clear. In [3^1 a UV 
cut off is imposed and then the glueball masses calculated - they are rather dependent on the 
choice of cut off. As noted by those authors the agreement with the Klebanov Strassler results 
is not so good. 



u 



III 1 


I 








1 60 
/ 40 
/ 20 






-3 


-2.5 


f 


-1 . 5 


-1 


-0.5 














-20 






III 


IV 











Figure l:The glueball potential in M = 1* with a = 1 and: I) b = 0.1; II) b = 0.19; III) b = 0.2; 

IV) b = 0.25. 
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Table 1: Spectrum of M = 1* glueball masses from supergravity with 6 = and b = 0.19, not 
taking into account the mixing. In the next column the M = 1 glueball spectrum obtained 
from Klebanov Strassler background is shown for comparison. The last column shows the 
analytic result obtained in taking into account mixing. In all cases the ++ mass has been 

scaled to 1. 



3.2 Gluino Bound States 

We now move on to investigate the masses of bound states of the fermions of the M = 1* theory. 
For simplicity we will set a = (b = 0) in what follows. To study bound states of the massless 
gluino we look at fluctuations of the scalar a. The corresponding potential ()12|) is 



, ,2m. 
2cosh(- ? =) 



(26) 



Numerically plotting this potential reveals it to be bounded giving a discrete spectrum with 
a mass gap. As in the case of the glueballs, we can find the eigenvalues of (fTT|) with the shooting 
method. The results are shown in Table 2 below, normalized to the lightest glueball mass above. 
The gravity dual describes a moduli space of vacua corresponding to different background values 
of a so we would expect to find a massless boson. In fact we do and the smallest eigenvalue is 
zero. 

Considerable work has been done by the DESY collaboration jTH] studying pure M = 1 
Yang Mills on the lattice. They reach the supersymmetric theory by tuning the gaugino mass 
to zero and to date their results for bound state masses are still somewhat away from the 
supersymmetric point. In the future though it should be possible to extract the supersymmetric 
values and it would be nice to have AdS predictions available for comparison. The lattice theory 
at finite iV has an anomalous U(l)# symmetry and a complete mass gap. On the other hand 
in the M = 1* theory that U(l)# symmetry is only an accidental symmetry in the IR in the 
weak coupling limit and is not present at all in the large N theory [3^ plus the gravity dual 
describes a massless bound state. The fermionic vacuum structure is therefore quite different 
and it is hard to compare the results for fermionic bound states. The best hope would be to 
compare the ++ glueball spectrum above. 

Finally it is also interesting to look at bound states of the massive fermions to see how 
decoupled these fields are from the dynamics. There is an SO (3) symmetry acting on the 
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massive fermions so it is sufficient to analyse bound states of a single massive flavour. We 
therefore subdivide the operator O m further introducing the scalars ii and v corresponding to 

o, = ^2, a = E^<- (27) 

i=3 

This enables us to study bound states of the massive fermion ip2 by looking at fluctuations 
of \i. In terms of these scalars the potential is [21] 

V = | [-5 + cosh(4/i) - 4 cosh(2//)] - cosh( v / 2z/) [cosh(2^) + 1] 

. (28) 

+^ [-3 + 2 cosh(2 v / 2z/) + cosh(4^) 

Of course, we are still interested in preserving Af = 1 supersymmetry in the background. If 
we set 



-m, 



m 

7z 



(29) 



then we recover the M = 1* potential with a = 0. The Schroedinger potential we are interested 
in for fluctuations in ll about this background is 



, ,2m, , ,4m. 
-2-2 cosh(^) + cosh(^ 
V3 V3 



(30) 



This is indeed a bounded potential. The mass spectrum is found by shooting and the results 
are displayed in Table 2. We find that the fermion bound states are a little heavier than the 
lightest glueball but far from completely decoupled. 

Since these states are not decoupled we must worry about mixing between them and the 
states whose masses we have calculated above. In [23 E2 122] the authors make progress in this 
direction by calculating the spectra through the analysis of the two point correlation function 
for the associated operator, using the holographic renormalization method [3U1 EH] E2] The 
glueball spectrum associated to the dilaton in the M = 1* background was found analytically 
in |2Sj solving the full equations of motion, that is taking into account the mixing of the dilaton 
with other scalars in the supergravity theory. The authors found a spectrum given by m 2 = 
(pL) 2 4(n+l)(n+4), with n = 0, 1, 2, . . .. This spectrum, indicated in Table 1 above, is consistent 
at the 5-10% level accuracy with our numerical result. Our results for the gaugino bound state 
agree perfectly with their analytic expression for the spectrum m 2 = (pL) 2 4(n + l)(n + 2) with 
n = 0, 1,2, .... It is remarkable that the physics of the model is not dramatically affected by 
the mixing, and supports a more naive analysis where these extended super-partner states are 
neglected. 

It is interesting to note that if we study fluctuations of the supergravity scalar m which 
corresponds to a composite operator of all three fermions ()16j) we obtain the potential 



Un 



3 .a 9 . /. n 
2-" + A {A ] 



+ e 



2.4 



-2 cosh( 



2m ^ 
73' 



cosh( 



4m, 

V5' 



(31) 
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This potential is unbounded. This presumably corresponds to an instability in the field theory 
for such a bound state to decay to the states we have seen above. 
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Table 2: The first five bound states of the massless gluino ipi and the massive fermion -02 in 

the Af = 1* theory. 

4 The Yang Mills* Geometry 

The Yang Mills* theory ^3] is A/" = 4 super Yang Mills with mass terms for all the adjoint 
matter fields leaving just non-supersymmetric gauge fields in the IR. At large N the massive 
fields can again not be considered totally decoupled. 

The Yang Mills* geometry is obtained by turning on a supergravity scalar A from the 10 of 
SO (6) dual to the operator 

4 

= X>tfk. (32) 
i=i 

The potential for the scalar can be obtained from the M = 1* solution by setting the two 
scalars equal 1 : 

V = ~ (l + cosh 2 A) . (33) 
In this case M 2 = — 3 in (UJ) and the ultra-violet solutions are 

X = Mz + ICz 3 . (34) 

The field theory mass term has dimension 3. Thus in what follows M. = corresponds to a 
solution with just bi-fermion vevs while K — corresponds to the purely massive case. Giving 
a mass to all four fermions breaks supersymmetry completely so we would expect the scalars to 
radiatively acquire masses so that the deep infra-red should be pure Yang Mills. The analysis 
of the lOd lift in [T3| supports this hypothesis. 

lr To be precise one must set m — -y/3/4A and a — ^/l/4A to maintain a canonically normalized kinetic term 
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4.1 Numerical Solutions 



We are not able to write down first order equations for the deformations as there is no su- 
persymmetry and no superpotential. Thus we are forced to solve the second order equations 
numerically. Numerical solutions of the flow equations for the scalars \(z),A(z) are displayed 
in Figure 2 for different asymptotic boundary conditions. The mass only flow is unique, lying 
on the border between two separate behaviour flows. If there is even a small condensate X(z) 
diverges very rapidly. It is necessary to fine tune the initial conditions to high precision in order 
to isolate the mass only flow. All these flows are singular including the mass only solution. 

We need a criteria for deciding which of these flows is the physical flow for the YM* the- 
ory. The simplest criteria we have found is the boundedness of the ++ glueball Schroedinger 
potential. Thus we again look at linearized dilaton fluctuations and evaluate the Schroedinger 
potential in (|24j) . We plot the potential for a variety of flows with both mass and condensate 
present in Figure 3. The presence of a condensate leads to an unstable geometry, and an un- 
bounded glueball potential. For the mass only solution, however, the glueball potential is a 
bounded well and the spectrum is calculable. We therefore conclude that this is the physical 
flow. 

In [T]\ it was proposed that a necessary condition for singular flows to be physical is that the 
supergravity potential, evaluated on the solution of the equations of motion, should be bounded 
above. This condition is the origin of the bound (|23|). Figure 4 shows the supergravity potential 
evaluated for different asymptotic boundary conditions. All the Yang Mills* flows satisfy this 
condition, including the unphysical flows with a condensate. It is also a necessary condition for 
a confining gauge theory that the glueball potential be a bounded well. Here this appears to be 
a stricter condition that successful distinguishes the physical from the unphysical flows. 
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Figure 2: Plots of A and A vs z in the YM* set up for a variety of UV initial conditions 
I) condensate only (A ~ — z 3 ); II) mass only (A ~ —z); III) mass and condensate 

(\~-(z + z 3 )). 
The AdS function A is also plotted. 
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Figure 3: The glueball potential showing: I) mass only; II) mass and condensate; III) 

condensate only. 
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Figure 4: The supergravity potential for a range of flows: I) mass only; II) mass and 
condensate; III) condensate only. All cases are bounded above. 
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4.2 IR Asymptotic Solutions 

We have also been able to analytically isolate the IR behaviour of the unique YM* flow that 
gives rise to a bounded glueball potential. In the r coordinates, the infra-red corresponds to 
r — > r , A(r) — > oo and A(r) — > — oo. The flow equations become 

A" + 4A'A' = -- A e 2 \ 6A' 2 = A' 2 + ^e 2A , (35) 

and have solutions 

A = logh — - — | J , c =\nr> A = ^ l og\r -r \, (36) 
\|r — r |y V 9 3 

where r marks the position of the singularity. This is clearly not a complete set of solutions, so 
could represent any of the solutions in Figure 2. To check whether this flow really corresponds 
to the solution with just a mass term, we can compute the glueball potential. This can easily 
done by changing variables dz = e A dr 



Hence 



3 

and the glueball potential U (z) in the infra-red is given by 

U(z) = lA" + 9 -(A'r 



e- A dr (37) 
z-z = 3(r-r ) 1/3 . 



A(*) = 21og|l-^|, (38) 



z ) 2 ' 



(39) 



We see that U(z) — > oo in the infra-red. Comparison with numerical study of the full second 
order equations implies that these solutions are the IR limit of the mass only Yang Mills* 
geometry. If there had been a condensate present, we would have found U(z) — > — oo. 

4.3 Glueballs 

As discussed above, the ++ glueball Schroedinger potential for the mass only case is a bounded 
well. We calculate the glueball mass spectrum by shooting, and the results are shown in Table 
3. As in the previous example, the lightest glueball state is not a prediction, but can be used 
to fix the scale Aq C d- As in the Af = 1* case we normalize the lowest state to 1. 

Again a possible source of inaccuracy in the nuerical results can derive from the mixing of 
the dilaton with other scalars. Even if there is no analytic proof for such mixing, we expect it 
to occur by analogy with the supersymmetric case. 
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An AdS approach to non-supersymmetric Yang Mills glueballs has already been presented 
in the literature. The AdS Schwarzchild black hole background describing M5 branes with 
a compact time dimension describes the high temperature limit of a strong coupled 5d field 
theory [3J Ej. The effect of the temperature is to give masses to all fields except the spatial 
components of the gauge field. Thus in the deep IR this theory is a 4d Yang Mills theory. 
This geometry was used to compute the glueball spectrum in [2j and we display the results 
in Table 2 for comparison (an alternative calculation of the same spectrum is given in [37]). 
The results are surprisingly close, differing only at the 10% level. We believe that taking into 
account the mixing would increase the masses, as happens in the supersymmetric case, making 
the agreement even more good and more remarkable. Our theory is 4d M = 4 Yang Mills in 
the UV whilst the AdS Schwarzchild dual describes a higher dimension theory in the UV. Since 
the UV is not decoupled from the strong interactions we might have expected at least order one 
differences. The closeness of the results hints at the relative stability of glueball masses across 
a wide range of theories and strengthens the case for trusting the AdS method of computation. 

It is also interesting to compare the results to lattice computations in pure Yang Mills. Data 
is only available for the lightest two ++ states. In Table 3 we show the data for N = 3 [38J and 
also the iV — > oo extrapolation of jSH]- The lattice errors are of order 10%. The AdS results 
actually match better to the N = 3 results but still agree with the N — > oo result at the 15% 
level. The data is very limited but again the agreement is surprising given that the UV of these 
theories are all so different. 



State 


Yang Mills* 


AdS-Schwarz 


Lattice (N = 3) 


Lattice (N = oo) 


0++ 


1.0 (input) 


1.0 (input) 


1.0 (input) 


1.0 (input) 


++* 


1.5 


1.6 


1.5 


1.9 


Q + + ** 


1.9 


2.1 






qH — h*** 


2.3 


2.7 






Q-| 1-**** 


2.7 









Table 3: Spectrum of glueball masses from the lattice, from Yang Mills* and Witten's 
AdS-Schwarzchild dual. Again, the lowest glueball mass has been scaled to one in all cases. 

4.4 Fermion Bound States 

To test the degree to which the massive fermions of the YM* theory have decoupled we will 
calculate the mass spectrum of bound states for them. We use the M = 1* potential in (fT7|) 
where the scalar a describes the bi-fermion operator of a single flavour of adjoint fermions. 
There is again an SO (4) symmetry on the massive fermions so it is sufficient to study one 
flavour. Setting the fields m, o to their YM* background values we then look at fluctuations of 
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a. The appropriate potential (J 12)) for this scalar in the Yang Mills* background is 

V* = \a" + 9 -(Ay-3e 2A . (40) 

This is again a bounded well potential and the spectrum of bound states is calculable. 
The results are shown in Table 4. The bound states have masses of order the glueball masses 
and so can not be considered truly decoupled. Again this suggests that the success of the 
AdS computations of glueball masses to match lattice data above indicates the stability of the 
glueball masses across a range of gauge theories rather than that the decoupling of extra fields 
is being achieved. 



State 


mass 


1 


0.9 


2 


1.3 


3 


1.7 


4 


2.1 


5 


2.6 



Table 4: First five bound states of a massive adjoint fermion in Yang Mills*. 



5 Wilson Loops 

The quark anti-quark interaction potential may be studied in AdS duals |H| by introducing a 
probe D3 brane into the geometry at some radius z max . Fundamental strings between the probe 
and the central stack of D3 branes would represent W bosons which transform in the (N, 1) of 
the SU(N) x £7(1) gauge group - we may equally think of these states as quarks since they are 
in the fundamental representation of SU(N). Thus a string attached to the probe with well 
separated ends play the role of a quark anti-quark pair with mass of order of the energy scale 
determined by z max . The action of the string corresponds to the interaction energy between 
the pair. To study such a configuration we necessarily require the lOd lifts of our deformed 
geometries since the string lives in lOd. The lift of the Yang Mills* geometry has been found 
in ^3] while for the M = 1* geometry only a partial lift exists 12J. For these reasons we will 
work mainly with the Yang Mills* solution. Since the warp factor of the M = 1* theory behaves 
similarly to that of Yang Mills* it is probable that the two theories have qualitatively the same 
Wilson loop behaviour. 
The lOd lift of YM* is 

ds 2 10 = (£♦£-) *d*i 4 + (£+£-)-^4 (41) 
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where ds\ 4 is given in (Q) and 



ds 5 = £_ cos a dQ + + £ + sin a + £ + £,_da . 
The functions £± are given by 

£± = c 2 ± s 2 cos 2a, c = coshA, s = sinhA. 
The dilaton is given by 



The two-form potential is given by 



A( 2 ) = iA + cos 3 a cos 9 + d9 + A 



A_ sin 3 a cos 9_d9_ A 



with A± = sinh2 A/£±. Finally the four-form potential lifts to 



'(4) 



F + *F, F = cfe A dx 1 A dx 2 A fix 3 A rfcu 



(42) 



(43) 



(44) 



(45) 



(46) 



where u{r) = e 4A ^ A'{r). 

We can calculate the Wilson loop behaviour by lying a string in an xn direction and letting 
it move in r and a. The action for the string in Einstein frame is given by 



S 



2vra' 



d 2 ae^ 2 VdetG. 



(47) 



With this choice of orientation, the fundamental string does not couple to the background B- 
field (145(1 which has non-zero components only on the internal manifold. In the static gauge 



(48) 



(a = x 1 = x, r = x°) and using the z coordinates (ds 2 = e 2A ^ z \dx 2 + dz 2 )) we obtain 



S=L I dxe 2A Ji 
2tt ' v 



\ 



\dx \ dx 



where T is the period resulting from the time integration. The resulting equations of motion 
for z and a are then 



d_ 

dx 

d 

dx 



-2A 



i + z 



Vl + z' 2 + e~ 2A a' 2 



J2A 



d_ 

dz . 
d 

da 



-?A 



J2A 



e+V 1 + z' 2 + e~ 2A a' 2 



+ z' 2 + e~ 2A a 2 



0. 



0. 



(49) 



(50) 



Vl + z' 2 + e- 2A a' 2 _ 

These can be solved numerically but there are an array of solutions. The new feature 
relative to pure AdS is that the strings can have a non-trivial a profile. The reason for this is 
that the fermion masses we introduced break the SO (6) symmetry to 5*0(3) x 50(3) so there 
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are potentially different strings connecting particles in different subgroups of the SO (6). To 
begin with lets concentrate on strings that are between two identical particles and hence have 
no a variation. There are such strings - in the a EoM the last potential term is given by 



da 



,2A 



+ z 2 + e~ 2A a' 2 



sinh A sin 2a 

7^ 



-2A 



+ z' 2 + e~ 2A a' 2 



(51) 



so vanishes if a = nn/2. The kinetic term also vanishes if a =constant so these are solutions of 
the equations of motion. 

As a simple case we study strings for which a = n/2 which implies £+ = 1. The z EoM then 
becomes 



d 

dx 



e 2A z 



VTT7 2 



d(e 



2A\ 



dz 



l + z' 2 = 0. 



(52) 



This is straightforward to solve numerically in the background A(z) appropriate to YM* (shown 
in Fig 2). 

Consider initial conditions where we start the string at z — —0.5, which is in the AdS like 
region. We then vary the derivative of z' and shoot off strings to more negative z corresponding 
to the interior of the deformed AdS space - we show numerical results of this type in Fig 5. 
Initially as z 1 increases the strings penetrate the geometry more and return to z — —0.5 at 
larger x indicating that they describe a quark anti-quark pair that are more widely separated. 
This is standard AdS behaviour. However when the strings begin to enter the deformed space 
the behaviour changes. At a critical value of z' the string, although still penetrating deeper into 
the space, returns to z — —0.5 at a shorter quark separation. Thus there is a maximum quark 
separation described by the Euler Lagrange equation solutions. 




Figure 5: Wilson loops in YM* showing how the depth the probe string penetrates into the 
deformed space depends on the quark separation (or equivalently the initial condition of z'). 

To get a better handle on this behaviour we can use the x independence of the Lagrangian 
which implies the Hamiltonian is conserved 

V^e 2A 



VI + z' 2 + a' 2 e~ 2A 



constant. 



(53) 
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Thus, for a constant, 



o4A 



£4 



- 1. 



(54) 



We can determine the integration constant c by evaluating the above equation at the point of 
maximum depth reach by the string, zq, where z 1 — and which corresponds to x = (because 
of the the symmetry of the problem) 



c 2 = ^o)e 4A(2o) . 



The quark anti-quark separation is then given by 

tax 

dz 



L 

2 



zo 



„2 



- 1 



(55) 



(56) 



The energy of the string is given by S/T so we find 

E = — / dz 

IXC Jzo 



+ 



(57) 



These equations are again straightforward to solve numerically in the YM* background. In 
Fig 6 we show plots of the quark anti-quark separation vs Zq 7 the energy of the string vs z 
and finally the energy of the string vs quark anti-quark separation. We see again the maximum 
separation but now also that the strings that penetrate into the interior of the deformed space 
are higher action. 




1.75 -1.5 -1.25 -1 -0.75 -0.5 
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0.314 
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Figure 6: YM* Wilson loop results at a = 7r/2; we plot the quark separation L vs the 
maximum depth z that the string reaches; the energy of the string vs the maximum depth; 

and finally the energy vs the quark separation. 
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This behaviour naively suggests string breaking by quark anti-quark pair production (an 
interpretation we put forward in the first version of this paper). We must be careful though 
since the Euler Lagrange equations only return local turning points of the solution (one of our 
solutions for fixed separation must therefore be a maximum) and it is possible they miss some 
run-away direction in the potential. To address this issue we should consider a set of strings 
linking quarks at fixed separation but which penetrate the space to different degrees. Ideally 
one should include the true solution of the Euler Lagrange equations in this set but we will 
learn all we need from a simpler set of configurations. We take just a sinusoidal shape for the 
string 

z(x) = — Asimrx/L (58) 

where we vary the amplitude A. 

In Fig. 7 we plot the action of these configurations as a function of A for L smaller than, 
larger than and close to the critical quark separation found above. 




L(z=0.5) = 0.50 L(z=0.5) = 0.70 




L(z=0.5) = 0.86 L(z=0.5) = 1.4 

Figure 7: The action of sinusoidal string configurations against amplitude (penetration into 
the space) for four different values of the quark separation L. 

It is now clear what the Euler Lagrange equation solutions correspond to. There is first 
a minimum of the action corresponding to an AdS-like string, then a maximum (the second 
solution above) and finally the action falls as the string moves into the singularity at z — 
—2.3. When the quarks are close the global minimum is the AdS-like solution but at a critical 
separation there is a transition and the preferred configuration is that which falls into the 
singularity. We have also looked at solutions at fixed a = and with varying a. Although there 
is a little more structure in the a variation the same instability exists as that just described. 

We conclude that we can not understand the behaviour of the Wilson loops at the super- 
gravity level for the IR. One could imagine that with an appropriate stringy resolution of the 
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singularity, perhaps in terms of some fuzzy brane configuration, the situation might be reme- 
died. It could be that the strings prefer to form a bound state with the interior branes and lie 
along the singularity (the mechanism suggested by Polchinski and Strassler jTHj) giving rise to 
a confining potential but this must remain speculation in these cases. 

6 Summary 

In this paper we have studied aspects of confinement in two deformed AdS geometries, Yang 
Mills* and M = 1*. We have shown that they describe discrete glueball spectra indicating that 
the geometries do confine. Study of bound states of the massive fermions in the geometries show 
that there is little decoupling of these massive states. Nevertheless the glueball spectra agree 
well with other gravity duals with the same infra-red degrees of freedom, even if the ultra-violet 
of the theories differ significantly. 

We also studied the behaviour of Wilson loops in the lOd lift of Yang Mills*. We found 
though that when a quark pair is separated beyond some critical separation the geometry 
prefers a string configuration connecting the quarks that falls into the singularity. It is therefore 
not possible to understand Wilson loop behaviour without a stringy resolution of the interior 
which currently is lacking. 
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